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We discuss the connection between absorption, averages and fluctuations in nuclear 
reactions. The fluctuations in the entrance channel result in the compound nucleus, 
Hauser-Feshbach, cross section, the fluctuations in the intermediate channels, result 
in modifications of multistep reaction cross sections, while the fluctuations in the 
final channel result in hybrid cross sections that can be used to describe incomplete 
fusion reactions. We discuss the latter in details and comment on the validity of the 
assumptions used in the develpoment of the Surrogate method. We also discuss the 
theory of multistep reactions with regards to intermediate state fluctuations and the 
energy dependence and non-locality of the intermediate channels optical potentials. 
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I. INTRODUCTION 



In the theory of nuclear reactions, one relies on statistical ideas to simplify their many- 
body nature. The first step is to introduce average amplitudes. The second step is to 
devise a Schrodinger equation that supplies these amplitues (the optical model), and 
the last step is the calculation of the fluctuation cross section. Several methods are 
available that supplies the derivation of this fluctuation cross section. The concept of 
energy averages and/or ensemble averages are used to obtain the Hauser-Feshbach cross 
section and corrections to it. The ergodic theorem is commonly employed to argue the 
equivalence of the two averaging methods. As much as fluctuations arise when one uses 
average amplitudes, the concept of fluctuations can be generalized to multistep processes, 
where the intermediate channel green functions are replaced by their averages. Similarly, 
if an exit doorway dominates the reaction, such as the excitation of a giant resonance 
in one of the nuclei participating in the reaction, or the coupling to the continuum, 
conveniently discretized, the final channel fluctuations set in resulting in an incomplete 
fusion cross section In both of these cases, the absorption and fluctuation analysis is made 
on the exact entrance channel and exit channel wave functions. 

In multistep reactions, another important reaction operator comes into the picture. 
The many-body intermediate channel Green's function. Again, one can decompose this 
quantity into an optical piece and a fluctuation piece. The fluctuation part can be 
calculated using the different optical quantitites as they appear in the KKM Optical 
Background Representation. One important addition to this procedure is the inclusion 
of other collective states which might be excitated in the intermediate propagation of 
the system. We have in mind the Brink- Axel mechanism. This has been accomplished 
recently in the theoretical description of the excitation of multiphonon states. It is 
important to extend this new feature to multistep direct reactions in general. 

In this talk I will discuss the incomplete fusion reaction theory and its connection to 
the Surrogate Method. I also discuss the Brink-Axel effect in multistep reaction theory. 
In both cases the Absorption- Fluctuation "theorem" is invoked to get a practical closed 
theory of the reactions. Some comments concerning the foundation of the Statistical 
Multistep Reaction Theory are also made. 

The wave function fluctuation can be written as 

^ ^optical ^fluctuation 

where ^fluctuation can be related to the optical quantities contained in \& = ^ optical ■ 
By construction the energy average of the latter wave function is 0. This then results in 
an average cross section containing an optical one ( calculated with DWBA or Coupled 
Channels Theory) and a fluctuation cross section calculated with optical transmission 
coefficients ( or matrices). If the wave function refers to the final channel, then the 
resulting cross section will contain a direct one plus a fluctuation one which refers to the 
formation of a compound nucleus in a subsystem ( if one has in mind a breakup process 
a(b +x) + A — > b + (x +A) — > b + c + C. The total formation cross section of the 
compound subsystem, a(b +x) + A — > b + (x +A), is the incomplete fusion cross section. 

When dealing with multistep processes, the intermediate channel Green's function, Gi 
can be decomposed into an optical plus a fluctuation pieces. The optical part contains 
complex energies. The fluctuation one is a multistep Green's function itself. It is this 
latter which would, if the conditions are met, contain the Brink- Axel effect. 



II. FINAL STATE FLUCTUATIONS 

We consider the following reaction Q, H @, Q , B H 03 

a(x + b) + A -y b+ ^ ( X + A )i (!) 

allstates 

and treat b as a spectator, namely it only suffers elastic scattering in the optical potential 
field supplied by the target. 

The cross-section for observing b is ( no CN formation in the xA subsystem ), 

j (direct) r. 
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where p{E b ) is the density of final b-states, p(E b ) = ^pk* ' anc ^ ^ e " ne g a tive-energy 
entrance channel wave function" of x, pi + \r x ) is, 



P: 



W(r.) = (l^foMafa - r«)*i + W.)) (3) 



where is the full three-body wave function of the b + x + A system. Further, the 
imaginary potential W^\r x , Ei + B a — E b ) takes into account the direct reactions of x 
with the target. 

The above equation for the direct inclusive cross section cross does not take into 
account the fluctuation in the xA subsystem owing to the CN formation. The inclusion 
of this contribution can be easily made by adding another piece to the imaginary part of 
the xA optical potential. This is easily seen if we write the cross section for observing b 

as 



l^tE = J^ p{Eb) £ < * (+)(rx ' ^A\V bx \4- ] %A > KE -E b - E c ) 

b b a c 

'.(-) 



X 



<C ) ^AlHx|* (+) (r x ,r b )$ A >. (4) 

The above equation can be reduced using standard procedures into the following com- 
pact form 

j j elastic— breakup » inelastic— breakup 

da b _ da b F + da^ _ 



dVt b dE b dVt b dE b dVt b dE b 

where the inelastic breakup cross section can be further decomposed into, 

^ inelastic— breakup ^ (direct) ^ {incomplete— fusion) 

dQ b dE b dQ b dE b + dfl b dE b ' ^ 

From the general form of the imaginary potential one can write the following for the 
incomplete fusion cross section which is the final state fluctuation contribution 

^^.(incomplete— fusion) ^ 



dVt b dE b ' V a 



p(E b ) V -Zof(Ei + B a - E b ) (7) 



C 



where a F c is the fusion cross section from channel c ( corresponding to an excited 
non-elastic channel of the xA subsystem), which is populated by a. 



o?=^<P™\W c \p™>. (8) 

with 

\p ( c +) >= G^%AP^ ] > • (9) 
Cleraly, Eq.(6) can be written more generally as, after introducing the xA fusion 
imaginary potential, W^. A , 

^^inelastic— breakup ^ 



dVthdEh hv, 



£-p(E b ) < pi +) (r,)|[^j;V) + wZ\r)]\pM(v x ) > (10) 



a 



Eq. (10), with W^ F ' associated with the annihilation component, was used to calculate 
the inclusive annihilation of very low energy antiprotons on deuterium [sj. Further, the 
above theory has been also used to treat incomplete fusion reactions involving weakly 
bound exotic nuclei [10]. Of course direct capture of x, without the concomitant excitation 
of the subsystem is the first term in the sum above. It is the above incomplete fusion 
cross section which is invariably referred to as hybrid, Trojan-Horse [ll[ or Surrogate 
1~2|. (l3| . cross section. In all of these applications only the direct capture is taken into 



account. Further, the x-particle is taken to be a geniune projectile, without due attention 
to the fact that its energy is dispersed by the internal, Fermi energy, inside a, the carrier 
of x. 



III. GENERAL CHARACTERISTICS OF HYBRID REACTIONS 

In the Surrogate method one uses a description of a "Desired reaction" , Dr in terms of 
a " Surrogate reaction" , Sr. The Sr is invariably of the type d+D — > b+B* — ► b+c+C+.... 
and the Dr of the type a + A^B*^c + C+ ....[14j. The compound reaction in the 
Dr is described by the Hauser-Feshbach theory, which supplies the following for the cross 
section, 

J,7T 

where a^ N is the partial fusion cross section ( CN formation cross section) and can 
be calculated from the optical model, while the function G^ N is the branching ration for 
the CN decay. This latter quantity is difficult to predict. For this, one relies on the Sr, 
which gives, 

Px(E) = E F ? N ( E > J -> *) " G C X N (E, J, tt), (12) 

where the function Ff N (E, J,n) is evaluated from direct reaction theory and is in- 
timately connected to the partial incomplete fusion cross section of Eqs.(7-9). What is 
usually done, however, is to assume ( J, 7r)-independent G CN and employ the Weisskopf- 
Ewing description of the Dr and the Sr (P X (E) = G^ N (E)). This prescription finally 
gives, 



<*{E) = < N {E)-P X {E). (13) 

To date, most of the effort in the application of the Sr relied on the WE approximation 
and Eq. (13). The Surrogate Method has been quite successful for a variety of Sr and 
it would be only natural to asky why, since several rather drastic assumptions (e.g. the 
WE approximation, and the use of on-shell description) are employed. In the following, 
I raise several theoretical issues. 

1- "The "incident projectile" energy, E a is dispersed by an amount directly related to 
its internal, Fermi, motion inside the real projectile d. The wave function of the captured 
particle a is given by the source wave function, pa (r ), which describes the incoming a 
off the energy shell. If the dispersion in the a energy is ignored as done in the Surrogate 
Method, one may wonder how large the off shell effects would be. 

2- The capture of a into the a — A compound nucleus B* proceeds in a multistep 
fashion. This is quite apparent from the general structure of the fusion cross sections 
appearing inside the sum above. 

3- The critical questions in Surrogate Methods are not the CN decay ratios. Rather, 
it is the formation of the CN. 

4- The d particle could be a " phonon-photon" as in (a, a') reactions. How to treat 
these in view of the theory above? One may env isag e using the equivalent nuclear-phonon 
method, as develpoed by Feshbach and Zabek [15| 



IV. INTERMEDIATE STATE FLUCTUATIONS 

In the excitation of nuclear states through multistep processes one may find it neces- 
sary to take into account the fluctuations inherent in the intermediate propagation of the 



system. This can be important in the excitation of multiphonon states [16l . 1 171 1 la . Il9l 
20| . If the collision time is comparable to the decay time of the intermediate collective 
state, then through the damping one may excite a collective Brink-Axel phonon on top 
of the background of the one phonon state. This will add a different type of fluctuation 
cross section, which has to be added to the usual mutistep cross section. As a result, 
the low-energy multistep direct cross section has to be modified. The basic quantity that 
averages and consequently develops fluctuations is the intermediate Green function. To 
see how the intermediate state fluctuations affects the cross section we first give a very 
brief description of the |2lJ statistical multistep direct reaction theory. 

The n-step transition from the initial channel i to the final channel / through the action 
of a transition potential operator V is given by an amplitude which has the general form, 

T if = V fn ^G { n + \(E) ■ ■ ■ V 32 G { 2 +) (E)V 21 G[ + \E)V U (14) 
where the intermediate channel Green's function G^(E) is given by, 

= E^TTe (15 > 

where Hj is the many-channel (body) Hamiltonian describing the colliding nuclear 
system in channel j. 

To be specific, we consider the excitation of multiphonon giant resonances. In this 
case the channels j refer to doorway channels, with d\ being the the first doorway or 



GDR, while d<i the second doorway or the Double Giant Dipole resonce (DGDR) etc. 

The propagator in the space of d\ is then denoted by the Green's function G d *(E). 
This Green function can be written as: 



G^{E) = G^{E)+G^ fluct \E) (16) 



where G { +\e) is the average Green function containing the damping width of the 
intermediate one-phonon collective state viz, 



1 



G%\E) = : — (17) 

E - H x - ex + i-£ 

and the fluctuation contribution contains explicit reference to the fine structure states 
1 

di 



that give rise to 



G^ fluct -(E) = G qi (E)v gidl G<£\E) (18) 

The Green function G qi represents propogation of the system in the qx subspace, 
the latter being spanned by the fine structure states to which the first doorway ( the 
single-phonon resonance) is damped and as a consequence it acquires the damping width. 
When calculating the two-step cross section, |0 >— > \dx >— > |c?2 >, one would obtain the 

average one, related to G%\e) pi us a fluctuation one related to G { d )fluct \E). This latter 
can be evaluated using the usual energy averaging procedure of products of two rapidly 
fluctuating Green functions; G qi (E)G q , (E f ). The resulting cross section is proportional 

to F d T c (E), where t c (E) is the collision time. If we call the average cross section of the 
two-step process, W^ 2 \ the usual multistep (MS) cross section for going from the ground 
state to the one-phonon doorway followed by the transition to the final, two-phonon 

doorway, through G d ^(E), by crjj/ s , and the two-step fluctuation cross section, by aff , 
then we find, 

(2) _ ^(E) (2) 

a fi — 2% MS ^ ' 

Accordingly the energy averaged cross section is 

wm= ( l+ <^Eiy aSs (20) 

It is evident that the fluctuation contribution to the two step cross section resulting 
from the excitation of the Brink- Axel phonon on top of the background of the one-phonon 
doorway, could become important at energies where the collision time is appreciable ( low 
energies). In the case of the excitation of the double giant dipole resonance (DGDR) in 
the reaction 208 Pb + 208 Pb at 640 MeV.A bombarding energy the fluctuation contributio 



is found to contribute 30% to the cross section |17l |. At the lower energy of 100 MeV.A 



the contribution is 100%. The same findings were reported in [20| where the method of 



supersymmetry ensemble average [23J was employed to obtain the Brink- Axel fluctuation 



term in the two-step cross scetion. The method of [18| can also be used to derive the 



cross section for the three-step (TGDR) process Ji| and the general form of the n-step 
cross section, viz 



a 



(n) 



n 



k)\ (n 
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(n + k) 




MS 



:2V. 



It is therefore clear that intermediate state fluctuations must be given due considerations 
in the theory of multistep direct reactions. It remains to be seen what has to be done in 
the SMSD theory of in order to take into account the intermediate state fluctuations 
and assess the convergence of the resulting modified expression of the cross section. 



V. NON-LOCALITY, ENERGY-DEPENDENCE AND COMPLEXITY OF 

THE OPTICAL POTENTIAL 

The optical potential is the backbone of nuclear reaction theory. The original mo- 
tivation behind the introduction of the OP was the need to analyse neutron scattering 
data with the one-body Schrodinger equation. The general properties of the OP can be 
summarized by three basic characteristics: Energy-dependence, Non-locality and Com- 
plexity. All of these features are clearly exhibited in the Feshbach form of the Optical 
Potential Operator, 

PU{E)P = PVQG QQ {E)QVP (22) 

where V couples the P-space to the Q-space and Gqq(E) is the full Green function in 
the Q-space. Here the P-space and the Q-space could be the open channels and closed 
channels (CN) subspaces, respectively, or simply two sets of open channels. In the case 
that Q represents closed channels, one resorts to energy averaging to smooth out the 
fluctuations in Gqq(E) due to the CN resonances. This can be expediently accomplished 
with the Kerman-Kawai-McVoy Optical Background Representation [221 ] . The energy 
averaged Gqq(E) is then just 

G^{E) = G QQ {E + iI) (23) 

where I is the averaging energy interval. It is then evident that the complexity of PU {E)P 
arises from the average intermediate state propagation in the Q-subspace. The optical 
potential is non-local owing to the same intermediate state propagation, 

PU(E, r, r')P = PV(t)QGqq(E + U)QV{r')P (24) 

If the Q-subspace represents other open channels not contained in P, the same arguments 
as above apply, with the Q-channel propogator now acquiring an imaginary part owing 
to the boundary condition of open channels ( the ie factor). 

Evidently, non-locality, complexity and energy-dependence of the Feshbach optical 
potential originate from the same source; the intermediate, Q-space, propagation of the 
system. Of course, the non-locality arising from the exchange effects should also be 
present, through the addition to the Feshbach potential, what is commonly called the 
bare potential, which is taken to be real. This is the mean field potential acting at positive 
energies. In actual applications in nuclear reactions of relevance to this Workshop, the 
optical potential is usually taken to be local, complex and with slow energy dependence. 



Obtaining an equivalent local potential from the non-local one is justified using the Perey- 
Buck|24] prescription, which brings in a non-despersive energy dependence besides the 
dispersive one contained in the intermediate Q-space propagator ( see above ). 

The presence of energy dependence in the intermediate channel optical potential which 
appears in the Green function brings in problems related to the contruction of the dual 
states in the spectral decomposition of the latter. These dual scattering states can not be 
obtained from the optical Schrodinger equation by merely changing the optical potential 
by its complex adjoint as is customary {25|, (26[, [23]. In fact, the dual scattering states for 
energy-dependent optical potentials can only be formally obtained by solving an integral 
equation containing in its kernel the half-on-energy-shell physical T-matrix. Of course 
when formulating the M SD reaction without the use of the spectral decomposition of the 



intermediate channel Green function as was done in [28|, [29| for two-step processes, the 
above question of the dual state is avoided. On the other hand, to extend such theories 
to processes involving more steps is rather impractical and the use of procedures such 
as those of [21] become unavoidable as long as the correct handling of the dual states 
is employed. First we recall the approximations employed by [2l[ to obtain, from Eq. 
(14), a cross section for multistep processes which is given by a convolution of single 
step, DWBA-like, cross sections. Let us suppose that the intermediate channel j Green's 

function is replaced by its energy averaged part, and call it G^ + \E). This Green's 
function is given by 

1 



In the above, the optical potential Uj(E) is, as already alluded above, energy- 
dependent and consequently complex, non-local potential. The three basic assumptions 
of [2l| are: a) ingnore the energy-dependence of the optical potential but keep it complex, 
b) use the spectral representation of the resulting Green's function in the form, 

(+) ,?,(+) 1 



and, c) take only the on-energy-shell part of the above ( the delta function part), viz 



(27) 

h 2 k 2 

where Ej = E — Ej = -^- L 

With the above form of the Green's function, the multistep series can be collapsed 
into a product of one-step amplitudes of the type, 

T 12 =<¥£\V 12 \^ > (28) 

The above amplitude is not of the usual DWBA type as the dual solution < | 

appears instead of the usual < ^ |- Oily the last term in the series involving the 
transition to the final channel is of the usual DWBA form. In the original FKK paper 
211 such a distinction was not made and all terms in the series were taken to be of the 



DWBA type. This assumption, however, was corrected later by [2 



In the following we argue that the energy dependence of the intermediate channel 
Green's function, even if it is weak, does not allow one to use the dual wave function 
above, obtained in [26[ and [27[ by solving the optical Schrodinger equation with Uj 
replaced by its complex adjoint llj. 

Consider the Schrodinger equation with energy-dependent interaction 

(E k -H -U(E k ))\¥+ ] >=0 (29) 
The dual scattering state is then defined such that 

< $L +) |^q +) >= (2vr) 3 5(k - q) (30) 

The equation that determines the dual state is then easily derived from the Lippmann- 
Schwinger equation of the physical state and the above orthonormality condition 

< ^ +) |q > + < ¥ k +) \Gi + \E q )U(E q )\^ >= (2nf6(k- q) (31) 

The above equation reduces to the usual one employed in [26[ in defining the mod- 
ified DWBA amplitudes in the FKK, MSD cross section, if U is taken to be energy 
independent, viz 

< ^ +) |q > + < ^ + VG«(^)|q >= (2n) 3 5(k- q) (32) 

In general, however, the energy dependence of U has to be taken into account and 
Eq. (21) has to be dealt with. To exhibit the physical quantities that appear in the 
integral equation that determins the dual sate, we use the spectral decompsotion of the 
free Green function and use the definition of the half-on-energy- shell T-matrix, 

<k'|r(^)|q>=<k'|C/(^)|*W>, (33) 

to obtain 



< ^ > + / < ¥ ^ > E ,-l + u < k 'i T < £ °>i q >= (^("-o) < 34 > 

It is thus clear that one needs the complete knowledge of the half-on-enegy-shell phys- 
ical T-matrix in order to solve the integral equation to obtain the dual state. Even if 
it is obtained, one would have to worry about the spectral decomposition of the inter- 
mediate channel Green function. For this reason, I would suggest the introduction of an 
equivalent energy-independent optical potential, U which, though highly non-local, can 
be employed in an integrodifferential equation for both the physical state ( with U) and 

the dual state ( with U ). The definition of U is just 

U(E k )\¥ k +) >=U\*l +) > (35) 

The obtention of U has been discussed during the 80's. I refer to [3(| for a dicussion of 
the contruction and limitation of such an equivalent energy-independent optical potential. 



VI. CONCLUSIONS 



Absorption and fluctuations go hand in hand in nuclear reactions. Initial state fluctu- 
ations give rise to the compound nucleus, Hauser-Feshbach, cross section. Intermediate 
state fluctuations could add contributions to multistep processes. Final state fluctua- 
tions give rise to incomplete fusion. In this talk I have elaborated on these effects. In 
particular, the Surrogate Method, which is incomplete fusion followed by the decay of 
the xA subsystem, is discussed and means to assess its limitations have been laid out. 
The question which theorists should address is why does this method work so well, as we 
have heard from several speakers in this Workshop. Further, the multistep direct reaction 
theory has also been considered. I have pointed out means to improve the convergence 
of the FKK version of the theory by including the intermediate state fluctuations and 
through a more careful treatment of the intermediate channels' optical potentials and the 
dual states. 
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